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1 Introduction 



Interactions of hadrons at high energy — such as those probed in ep scattering at HERA, pp 
scattering at the Tevatron, and pp scattering at the forthcoming LHC — are to be understood 
in terms of the interactions of their quark and gluon constituents. The parton distribution 
functions (PDFs) that describe the proton's quark and gluon content are therefore essential 
for testing the Standard Model and searching for New Physics. 

The PDFs are functions f a (x,fi) where x is the fraction of the proton's momentum 
carried by parton species a at scale /x, in a frame where that momentum is large. For small 
values of /x, corresponding to long distance scales, the PDFs express nonperturbative physics 
that is beyond the scope of present calculations from first principles in QCD (although some 
progress has been made using lattice methods [1]). Current practice is instead to parametrize 
the PDFs at a scale /x that is large enough for f a (x, /x) to be calculated from f a (x, /x ) at all 
\i > fiQ by perturbation theory. The unknown functions f a {x, /Jo) are determined empirically 
by adjusting their parameters to fit a large variety of data at /x > jiQ in a "QCD global 
analysis" [2,3]. 

A number of important processes, including Higgs boson production in certain scenarios, 
are particularly sensitive to the bottom or charm quark distributions f c (x, /x) and fb(x, /x). In 
the global analyses that have been carried out so far, it is assumed that the charm content of 
the proton is negligible at /x~m c , and similarly that bottom is negligible at /x^mj, so these 
heavy quark components arise only perturbatively through gluon splitting in the DGLAP 
evolution. The global fits are not inconsistent with this assumption, but the data sets they 
are based on do not yet include experiments that are strongly sensitive to heavy quarks, 
so substantially larger c or b content cannot be ruled out. Direct measurements of c and 
b production in deep inelastic scattering are also consistent with an entirely perturbative 
origin for heavy quark flavors [4], but those experiments are not sensitive to heavy quarks 
at large x. 

Meanwhile, in the light-cone Fock space picture [5], it is natural to expect nonpertur- 
bative "intrinsic" heavy quark components in the proton wave function [6, 7]. Furthermore, 
s and s quarks each carry ~1% of the proton momentum at /xo = 1-3 GeV [2], which implies 
that states made of uudss, together with gluons, make up a significant component of the 
proton wave function. By analogy, one would expect uudcc and uudbb also to be present— 
although the degree of suppression caused by greater off-shell distances is difficult to predict. 
A suppression as mild as ~ has been derived using a semiclassical approximation for 

the heavy quark fields [8]. 

An alternative way to describe the proton in light-cone Fock space is in terms of off-shell 
physical particles — the "meson cloud" picture [9-11]. Specifically, the two-body state D A+, 
where D is a uc meson and A+ is a udc baryon, forms a natural low-mass component. This 
is the flavor SU(A) analog of the K + A° component that is a natural source of strangeness, 
and in particular of f s (x,fi) ^ f s (x,fi) [12]. A charm contribution from the two-body state 
p J/0 is also possible. 
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The light-cone view is not developed to a point where the normalization of uudcc and 
uudbb components can be calculated with any confidence, though estimates on the order of 
1% have been found for intrinsic charm (IC) using a meson cloud model [9], the MIT bag 
model [13], and an SU(4) quark model [14]. However, we can use the picture to predict 
the x-dependence of the non-perturbative contribution. A central feature of the light-cone 
models is that heavy quarks appear mainly at large x, because their contribution to the off- 
shell distance is proportional to (p]_ + m 2 )/x, so the suppression of far off-shell configurations 
favors large x when m is large. We will show that this feature leads to similar predictions 
for the shape in x from a wide variety of specific models. 

Using the rough consensus of the models as a guide to the shape of x-dependence 
for intrinsic charm and bottom, it will be possible to estimate their normalization from a 
limited set of data. This will be carried out in a future publication. When more complete 
data become available, such as jet measurements with c- and 6-tagging (either inclusive jets 
or jets produced in association with W, Z, or 7), it will become possible to extract the 
x-dependence empirically. It will then be interesting to see if the model predictions for the 
x-dependence are borne out. 

Models in which the uudcc Fock space component is considered directly are described 
in Sec. 2. Models based on low-mass meson+baryon pairs are described in Sec. 3. The 
model results are compared with expectations for perturbatively generated heavy quarks in 
Sec. 4. The model results are compared with the light quark and gluon distributions in 
Sec. 5. Conclusions are summarized in Sec. 6, where simple parametrizations of all of the 
model predictions are tabulated for convenience in later work. The connection between the 
light-cone description and ordinary Feynman diagrams, which is used in Sees. 2 and 3, is 
derived in an Appendix. 

2 Five-quark models 

The probability distribution for the 5-quark state uudcc in the light-cone description of the 
proton can be written as 

d p=M ») iW* (2) d>x) J7 1 ^ , (i) 

j=1 Xj j=1 j=1 j=1 I s m o) 

where 

5 

s = 52(pL + m f)/ x J ( 2 ) 

and M is a normalization constant. Eq. (1) contains a wave function factor F 2 that charac- 
terizes the dynamics of the bound state. This factor must suppress contributions from large 
values of s to make the integrated probability converge. An elementary derivation of Eq. (1) 
is given in the Appendix. 
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2.1 The BHPS model 



A simple model for the x-dependence of charm can be obtained by neglecting the p± content, 
the 1/xj factors, and F 2 in Eq. (1). Further approximating the charm quark mass as large 
compared to all the other masses yields 



dP a II dx i S ^ ~ x ^ ( l / Xi + l / x ^ 

3=1 3=1 



(3) 



where £4 = x c and £5 = x s - Carrying out all but one of the integrals and normalizing to an 
assumed total probability of 1% yields 



dP 

dx 



fc(x) = fc{x) = Qx 2 [Qx(l + x) lnx + (l-x)(l + 10x + 



X 



(4) 



where x = x c or x 5 ■ Equation (4) was first derived by Brodsky, Hoyer, Peterson and Sakai 
[6], and has been used many times since. We will use this BHPS model as a convenient 
reference for comparing all other models. 

Charm distributions that arise when the transverse momentum content of Eq. (1) is not 
deleted are derived in the following subsections. 
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Figure 1: Momentum distribution of c or c from the 5-quark model with the exponential 
suppression of Eq. (5) (left), or the power-law suppression of Eq. (6) (right). Solid curves 
are the BHPS model of Eq. (4). All curves are normalized to 1% integrated probability. 



2.2 Exponential suppression 

A plausible conjecture would be that high- mass configurations in Eq. (1) are suppressed by 
a factor 

F 2 = e" (s_m o)/ A2 . (5) 



This exponential form makes the total probability integral converge for any number of con- 
stituents, while a power law would not (see Eq. (21) in the Appendix). 

Figure 1(a) shows the charm distribution for several choices of the parameter A in 
Eq. (5). The mass values used were mo = 0.938 GeV and m± = m 5 = 1.5 GeV. Constituent 
quark masses ni\ = m2 = = mo/3 were used for the light quarks, but even setting those 
masses to zero instead yields very similar results. All curves are normalized to 1% integrated 
probability. The results are qualitatively similar to the BHPS model, but are somewhat 
smaller in the region x > 0.5, which is the most important region as shown in Sec. 5. 

2.3 Power-Law suppression 

Alternatively, we might assume that high-mass five-quark states are suppressed only by a 
power law, say 

F 2 oc (s + A 2 )"™ . (6) 

Figure 1(b) shows the charm distribution for several choices of the parameter A with n = 4. 
The results are again rather similar to the BHPS model, and again smaller than that model 
at large x. This behavior is fairly insensitive to the choice of n: similar results were found 
for n = 3 (not shown), while large values of n revert to the exponential form of Section 2.2. 
Values n < 2 are unphysical, since they would make the total probability diverge. The value 
n = 3 is perhaps the most natural, since it leads to a dependence ~ 1/m 2 that is in line with 
the result of [8]. 



Quasi-Two-Body 
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Figure 2: Momentum distribution of c and c from the (udc)(uc) model (7), and c = c from 
the (uud)(cc) model (10). 
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2.4 Quasi-Two-body suppression 



Another approach to the suppression of high-mass Fock space components can be made 
on the basis of quasi-two-body states, such as those that will be considered explicitly in 
Section 3. For instance, we might assume the relevant 5-quark configurations are grouped 
as (udc)(uc), in which case a plausible wave function factor would be 

F 2 oc (s 124 + A^)" 2 x (s 35 + A 3 2 5 )" 2 (7) 

where 

S124 = (j>l± + m l)/ X l + (pL + m 2)/ X 2 + (P4x + m 4)/ X 4 (8) 

S35 = (pL + m l)l x z + (pL + m l )/ X 5 ■ (9) 

Figure 2 shows the c and c distributions according to this assumption. The parameter choices 
were A 12 4 = 2.5 GeV and A 35 = 2.0 GeV, but other plausible choices lead to similar results. 
Note that there is a small difference between the c and c distributions in this model, with — 
perhaps surprisingly — c(x) > c{x) at large x. 1 Observing a c(x)-c(x) difference would of 
course definitively prove a non-perturbative origin for charm. 

Alternatively, we might assume the relevant 5-quark configurations are grouped as 
(uud)(cc), so a plausible wave function factor would be 

F 2 oc (s 123 + A 2 3 )- 2 x ( S45 + A 2 5 )" 2 . (10) 

The result of this assumption with A123 = 1 GeV and A45 = 3 GeV is also shown in Fig. 2. 
It happens to be very similar to the average of c and c from the preceding model. 

3 Meson+Baryon models 

Another way to picture the proton in light-cone Fock space is as a superposition of config- 
urations of off-shell physical particles. In particular, the two-body state D A+, where D is 
a uc meson and A^~ is a udc baryon, is a natural low- mass component. 

We can model the D A+ probability distribution in the proton using Eq. (21) with 
N = 2 and F 2 oc (s^a + A 2 ) -2 . The physical masses are m = 0.9383, mi = 1.8641, and 
m2 = 2.2849 in GeV. We then model the udc distribution in A^" similarly, using = 3 and 
F 2 oc (s udc + A|)~ 2 , with m = 2.2849, mi = m 2 = 0.938/3, and m 3 = 1.6. We similarly 
model the uc distribution in D using N = 2 and F 2 oc (s u5 + A^) -2 , with mo = 1.8641 and 
m i — m 2 — 1-60 . (The charm quark mass here must be taken > 1.55 to keep m 5 + m u > m-jy 
for stability.) 

1 One might have expected c(x) > c(x) because the c quark comes from the heavier (baryonic) subgroup. 
But in fact, the two subgroups share the proton momentum fairly equally, while the c retains more of the 
momentum of its (mesonic) subgroup because it shares that subgroup momentum only with a single quark. 
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Figure 3: (a) Momentum distribution of A+ from p — >• D A+ with A p = 10 (dashed), 4 
(solid), 2 (dotted), (b) Momentum distribution of c from A c — > udc with A a = 4 (dashed), 
2 (solid), 1 (dotted), (c) Resulting distribution of c in p with A p = 4 and Aa = 2 and BHPS 
model (solid). 



The c and c distributions in the proton follow from convolutions of the distributions 
defined above: 

dx 1 f 1 (x 1 ) / dx 2 f 2 (x 2 )8(x-xix 2 ) = / — fi(y)f 2 (x/y) . (11) 
Jo Jx V 

Figure 3 shows the c distribution from p — > -D°A+. The result is very similar to the BHPS 
model. Figure 4 shows the c distribution from p — > D A+. On comparing Figs. 3 and 4, we 
once again find c(x) > c(x) at large x — for the same reason as described in footnote 1. This 
was observed previously in a meson cloud model that is rather similar to this one [11]. It is 
opposite to the asymmetry predicted by [15]. 

A contribution from the two-body state p J/ip is also possible. It is even slightly favored 
over D Ajt by having a lower threshold mass: m p + rrtj/^ = 4.035 GeV < mo + "iA c = 
4.149 GeV. (This is in contrast to the SU(4)-analog case of strangeness, where K + A° is 
strongly favored over p0° by m K + + m A o = 1.609 GeV m p + = 1.958 GeV.) Figure 5 
shows the c = c distribution from the model of a p J/ip Fock space component and the model 
of c or c in J/0. A range of reasonable parameters for the suppression of large masses was 
tried and all give similar results. 




4 Comparison with perturbative cc and bb 

When normalized to 1% probability, the BHPS model predicts that a fraction J [fc( x ) + 
f 5 (x)]xdx = 0.0057 of the proton momentum is carried by nonperturbative charm. The 
models of Sections 2-3 give quite similar values, ranging from 0.0046 to 0.0073 . 

These possible intrinsic momentum fractions can be compared with the standard per- 
turbative contributions to the proton momentum, which are shown in Fig. 6 as a function 
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Figure 4: (a) Momentum distribution of D° from p — > D°A+ with A p = 10 (dashed), 4 
(solid), 2 (dotted), (b) Momentum distribution of c from D — > cu with A^> = 4 (dashed), 
2 (solid), 1 (dotted), (c) Resulting distribution of c in p with A p = 4 and = 2. 

of [i. (These were calculated from the CTEQ6.1 global analysis [2], with uncertainty ranges 
based on the eigenvector uncertainty sets [2, 16].) Note that the c + c and b + b fractions have 
been multiplied by 10 for clarity. We see that a possible 1% intrinsic charm contribution 
would be rapidly overtaken by perturbatively generated charm, once the evolution in [i has 
proceeded a short distance above m c . 2 Gluon splitting similarly generates bb pairs rapidly 
above the scale yU~m{,. Thus intrinsic c and b cannot be expected to add significantly to the 
perturbatively generated c and b for most regions of x and /x. 

Nevertheless, the intrinsic cc component may be very significant at large x. This is 
demonstrated by Fig. 7, which shows the probability distributions as a function of x, weighted 
by a factor x 2 to clarify the large- x region. The intrinsic component is stronger than the 
perturbative one at x > 0.3, even for /i as large as 100 GeV. (The intrinsic component will 
of course also evolve with fi, but that will not significantly alter this comparison.) 

5 Comparison with light quarks and gluon 

Figure 8 compares the BHPS model, which is representative of all the models described here, 
with the light-quark flavors and gluon from CTEQ6.1. It shows the remarkable result that 
with the assumption of 1% intrinsic charm, the c and c distributions are larger than d and u 
at large x. This result arises from the m 2 /x term in the off-shell distance of the heavy-quark 
states. 

Figure 8 also shows that intrinsic c and c are much smaller than the valence quark 
and gluon distributions. This implies that the possibility of IC does not significantly affect 
the light quark and gluon distributions. As a corollary, it negates a pretty speculation 

2 The rapid rise of c + c is likely to be somewhat exaggerated in this CTEQ analysis, which uses the 
standard "zero mass scheme" wherein the charm quark is treated as a massless parton at scales fi > m c . 
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Figure 5: (a) Momentum distribution of J/ip from p — > p J/ip with A p j/^ = 5 (dashed), 3 
(solid), 1 (dotted), (b) Momentum distribution of c or c from J/ip — > cc with A gc = 5 
(dashed), 2 (solid), 1 (dotted), (c) Resulting distribution of c or c in p with A p j/^ = 3, 



that IC could be the source of an unexpected feature of the CTEQ6.1 PDFs which can be 
seen in Fig. 8: the gluon distribution is larger than the valence quarks at very large x for 
small /i. (That feature is not a robust feature of the CTEQ PDFs, since it can be made 
to disappear by a small change in the gluon parameterization at po, with an insignificant 
increase (A% 2 « 5) in the x 2 ~ 2000 of the global fit. The feature therefore does not actually 
require an explanation.) 

6 Conclusion 

The light cone ideas used here are at best qualitative and heuristic. It is not for example 
clear whether they should be applied in MS or some other scheme; or at what scale /io- They 
should nevertheless be a useful guide in the effort to measure intrinsic heavy flavors in the 
proton. 

We have shown that intrinsic charm (IC) will provide the dominant contribution to c 
and c at large x, if the shape of the IC distribution is given by the BHPS model and the 
normalization is anywhere near the estimated 1% probability. All of the other light-cone 
based models we have examined have roughly the same shape in x dependence, and hence 
they reinforce this result. Several of the models predict a difference between c and c, with 
c(x) — c(x) > at large x. Similar conclusions for the shape apply for intrinsic b. 

Assuming the 1% probability is approximately correct for IC, c is much smaller than u, 
d, and g at all x, so it has no appreciable impact on the evolution of other flavors. Intrinsic 
b is presumably even smaller. 

An estimate of (0.86 ±0.60)% was obtained for the IC probability some time ago [17, 18] 
by re-analysis of F 2 C data in deep inelastic muon scattering on iron [19]. That estimate 
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ju. [GeV] 

Figure 6: Fraction of proton momentum carried by valence quarks, gluon, b + b, and c + c 
as a function of scale /i, with uncertainty bands. The b + b and c + c curves were multiplied 
by 10 for clarity. The dot-dash line indicates the level of intrinsic c + c predicted by BHPS, 
also multiplied by 10. 

continues to be cited (see e.g. [5]) as evidence for the existence of IC, although it is obviously 
of limited statistical power; and when possible variations in the parton distributions are taken 
into account, the F 2 C data are consistent with no IC [11,21]. Measurements of _F 2 C at HERA 
[4] are also consistent with no IC, but those measurements are not at sufficiently large x to 
have any sensitivity to it. For other experimental indications of IC, see [20]. 

In order to actually measure intrinsic charm or bottom, it will be necessary to have data 
that are directly sensitive to the large-x component. Likely candidates are jet production 
with c- or 6-tagging — either inclusively or in association with W, Z, or high-p^ 7- It may 
also be possible to extract useful information from coherent diffractive dissociation processes 
such as p — > p J/ip on a nuclear target [22] . 

For convenience in future work, the model curves (1)— (12) that appear in Figs. 1-5 can 
be adequately represented by a simple parametrization which is given in Table l. 3 To create 
this table, the normalization coefficients A Q were chosen to make the momentum fraction 
Jo fc( x ) x or f f 5 (x) x dx equal to 0.002857, the value given by the BHPS model when that 
model is normalized to 1% probability. This is different from the normalization of the curves 
shown in the Figs. 1-5, which was such that each curve corresponded to 1% probability. 
This new normalization is more useful for applications, since it places more emphasis on 
large x, where intrinsic charm will be important if it is important at all. For comparison, 
the momentum carried by s or s at \x = 1.4 GeV is — as it should be — substantially larger 
(by a factor of 4) than this working estimate of 0.002857 for c or c. 

3 With the parameters listed for it, this simple parametrization also works very well for the BHPS model, 
though of course the full expression (4) is not inconvenient. 
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Figure 7: Distribution x 2 f 5 {x) from CTEQ6.1 (extrinsic charm) at fi = 2, 5, 100 GeV (short, 
medium, long dash) compared to BHPS model (intrinsic charm) (solid). The intrinsic charm 
component dominates at large x. 
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Appendix: Light-cone probability distributions from Feyn- 
man rules 

This Appendix shows how to derive light-cone probability distributions directly from Feyn- 
man diagram rules by a thought experiment. 4 

For simplicity, consider a spin particle with mass mo that couples to spin particles 
with masses mj, . . . by a point coupling ig, as illustrated in Fig. 9(a). The thought 
experiment consists of scattering this system at very high energy from a target that inter- 
acts with only one of the constituents, as illustrated in Fig. 9(b). This target supplies an 
infinitesimal momentum transfer that puts the iV-particle system on mass-shell (diffractive 
dissociation) with a cross section that must be proportional to the probability for that system 
in the original Fock space. 

Assume that the target provides a constant total cross section <jq, with transverse mo- 
mentum transfer dependence proportional to exp(— /3q±). The elastic differential cross sec- 
tion is 



da a, 



L e -Ptt (12) 



dq[ 167T 



The technique of this thought experiment can in fact be used to calculate coherent production such as 
occurs in the " A 2 / 3 " component of J/tfj production on nuclei - see [22] . 
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Figure 8: Parton distributions for u (long dash), d (short dash), u (long dash dot), d (short 
dash dot), and g (dotted) at fi = 2 GeV, compared with c from the BHPS model (solid). 




Figure 9: (a) Point coupling model; and (b) diffractive dissociation process used to compute 
its light-cone probability distribution. 



and hence the integrated elastic cross section is 



" = TSjT (13) 



2 



By Feynman rules, Fig. 9(b) gives an amplitude 



M= 2i9P*rQ { _ p £ /2) (14) 



and a cross section 



da = — d [ q M ff ( d3pj m , ) 5^(p + k-Y Pj -q)\M\ 2 . (15) 

Assume the gaussian parameter (3 is large, corresponding to a large spatial extent of the 
target in impact parameter. We can then set q± = everywhere except in the exponential 
factor and carry out the integral over q±. The Fock space probability density dP can be 
identified from the obvious relation 

da = cr cl x dP . (16) 
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Now introduce the light-cone components of all four-momenta, = ± p^) / a/2 , 



and define the light-cone momentum fractions Xj = Pj/Po- The components p^ and 
q(~^ are taken to be large, with p ± = q± = 0. The small components are determined by 
mass-shell conditions, e.g., p 2 = m 2 =>• = (p 2 ± + m 2 )/{2p^). This leads to 

^-(le^n(^) ^|>J*(i-|>)(^. (it) 

The covariant off-shell distance can be expressed conveniently in the form 

(TV-1 \ 2 
Po - 5^ I = Xn ( s ~ m o ) ' ( 18 ) 

where 

TV 

This leads to the final result 

TV TV , TV TV 2 

for the point-coupling model. Note that this result is completely symmetric in the particles 
1, . . . , N as it should be — it does not depend on which particle was singled out to scatter in 
the thought experiment used to derive it. It is straightforward to include more complicated 
vertices and factors due to spin using this Feynman diagram approach. When that is done, 
the result can depend on which particle is assumed to scatter, but the ambiguity vanishes 
at the pole at s — m 2 . Unitarity effects that keep the total probability equal to 1 could also 
be included. 

In this simple point-coupling model, high-mass Fock states are suppressed only by the 
"old-fashioned perturbation theory energy denominator" factor (s — ttlq)" 2 . To make the 
model more realistic, there must be a further suppression of high mass states associated 
with wave function effects — if only to make the integrated probability finite. It is natural to 
suppose that the additional suppression is a function of s. 5 

When a wave function factor [F(s)] 2 is included in Eq. (20), the transverse momentum 
integrals can be carried out by inserting the identity 1 = J S(J2ji.Pj± + m j)/ x j ~ s ) ds, and 
then Fourier transforming this delta function and the transverse momentum conserving one. 
The result is 



dP = I fr^U V,) r -*>)"-* ™? *> 



r- 



It is not natural to assume that the suppression is a function of the covariant variable t, as is done in 
some "meson cloud" models [23], since that assumption spoils the independence on which particle is taken 
to be off shell. A related argument leading to this conclusion is given in more recent meson cloud work [24] . 
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where sq = Ylj=i m j 'l x i • Note that [-F(s)] 2 must go to zero faster than l/s N 3 as s — > oo 
to make the integrated probability converge. 
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Model 


A 


Ai 


A 2 


(1) 


520.517 


4.611 


11.477 


(2) 


0.904 


1.271 


5.703 


(3) 


0.179 


0.496 


4.164 


(4) 


0.785 


1.220 


5.519 


(5) 


0.401 


0.886 


4.908 


(6) 


0.187 


0.521 


4.194 


(7) 


0.387 


0.822 


5.017 


(8) 


0.374 


1.010 


4.422 


(9) 


0.473 


0.994 


4.986 


(10) 


2.238 


1.897 


6.095 


(11) 


2.245 


2.511 


4.929 


(12) 


591.400 


5.065 


10.708 


BHPS 


1.052 


1.524 


5.377 



Table 1: Coefficients for a simple parametrization of the models: f c {%) — dP/dx = 
Aqx Ai (1 — x) A2 . The normalization Aq is chosen to make the momentum fraction 
Jo fc(%) xdx or f Q f 5 (x) xdx equal to 1/350, the value it has in the BHPS model. 
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